Painleve singularity structure analysis of three component 
Gross-Pitaevskii type equations 

T. Kanna* and K. Sakkaravarthi 
Department of Physics, Bishop Heber College, Tiruchirapalli-620 017, India 

C. Senthil Kumar ' 

Department of Physics, VMKV Engineering College, 
Periaseeragapadi, Salem-636 308, India 

M. Lakshmanan * 

Centre for Nonlinear Dynamics, Bharathidasan University, Tiruchirapalli-620 024, India 

M. Wadati § 
Department of Physics, Faculty of Science, 
Tokyo University of Science, 1-3 Kagurazaka, 
Shinjuku-ku, Tokyo 162-8601, Japan 
(Dated: October 26, 2009) 

Abstract 

In this paper, we have studied the integrability nature of a system of three coupled Gross- 
Pitaevskii type nonlinear evolution equations arising in the context of spinor Bose-Einstein con- 
densates by applying the Painleve singularity structure analysis. We show that only for two sets 
of parametric choices, corresponding to the known integrable cases, the system passes the Painleve 
test. 
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I. INTRODUCTION 



Integrable multicomponent nonlinear Schrodinger type equations have attracted con- 
siderable current interest in soliton research. Much focus has been paid to identify new 
integrable multicomponent type equations due to their many faceted applications in 
different fields of science such as nonlinear optics, Bose-Einstein condensates, biophysics, 
plasma physics, etc. 1-5 . Painleve singularity structure analysis is one of the powerful 
tools to isolate and identify integrable dynamical systems 6-10 . This procedure nicely 
complements other integrability tools like inverse scattering transform (1ST), infinite 
number of involutive integrals of motion, symmetries, Backlund transformations, Hirota's 
bilinearization method, etc., to study the integrability properties of nonlinear systems 1 ' 2 . 
By applying the Painleve test for integrability a class of integrable coupled nonlinear 
Schrodinger (CNLS) type equations, which arise in different physical contexts, has been 
identified 11-16 . 

In this connection, the system of CNLS equations in the presence of confining potential 
becomes the coupled Gross-Pitaevskii (GP) equations, governing the dynamics of two 
component Bose-Einstein condensates 17-19 . This kind of multicomponent condensates can 
also be created with the mixture of two different atomic species or by considering the 
hyperfine spin of atoms in the presence of optical dipole traps 20-22 . The latter entities are 
the so-called spinor Bose-Einstein condensates (BECs). 

Spinor Bose-Einstein condensates of ultra cold atoms can be created by liberating the 
hyperfine states by means of optical trapping. Two component condensates have been 
realized in 87 Rb (see Ref. 23) and also optically trapped three component condensates were 
studied in Refs. 24-27. The evolution of the spinor condensate wave functions is governed 
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by the following set of three-coupled nonlinear Schrodinger type equations , 

h 2 

ihi/j +ltT = -— ip+i,xx + (co + c 2 )(|^+i| 2 + \ipo\ 2 )i>+i 
2m 

+(c - c 2 )\i)-i\ 2 i)+i + c 2 v>*i^o> ( la ) 
h 2 

ih^ ,T = -7T- IpOXX + (Co + C 2 )(|^+l| 2 + |^-i| 2 )^ 

2m 

+c \ip \ 2 ^ + 2c 2 ^o^+i^-i 5 ( lb ) 

ihlp-i >T = -7T- fp-l,XX + (CO + C 2 )(|^-l| 2 + |-0o| )'0-l 

+ (C - C 2 ) | 2 V-1 + 02^+1^0 » ( lc ) 

where V±i,o' 8 are the wave functions of the three spin components, T is the time and X 
denotes the spatial co-ordinate. The effective one-dimensional coupling constants Cq and 
c 2 representing the mean field and spin exchange interactions, respectively, are given by 

Co = floi^a, c 2 = S^ffi where g f = ^ (i=c^) ' ^ = °' 2 " Here a /' s are the s ~ wave 
scattering lengths in the total hyperfine spin channel /, a± is the size of the transverse 
ground state, m is the atomic mass and the constant C = —£(1/2) ~ 1.46, where ( is 
the Reimann zeta-function. With the redefinition of T = ht, X — —&=x and transforming 
(ipi,ipo,ij)-i) — > ("01, v2-0oj'0-i)3 we can rewrite Eq. (CQ) in the standard form as 

iV>+i,t = -0+1,^ + (co + c 2 )(|^+i| 2 + 2|^ | 2 )^ + i 

+ (co-c 2 )|^_i| 2 V' + i + 2c 2 d^o, (2a) 

#0,i = -1p0,xx + (Co + C 2 )(|^+i| 2 + |^~l| 2 )^0 

+2co|^o| 2 ^o + 2c 2 V>oV>+i^-i, (2b) 
= -i>-i,xx + (co + c 2 )(|^-i| 2 + 2|^ | 2 )^-i 

+ (co-c 2 )|^+i| 2 V'-i + 2c 2 ^; i ^ 2 . (2c) 

We refer to Eq. (2) as the three-component GP type equations. The above system of 
equations has been solved by the 1ST method and multicomponent bright and dark solitons 
have been reported for specific choices of Co and c 2 29-31 . Now it is of interest to isolate all 
the possible integrable models arising from Eq. ([2]) for arbitrary choices of Cq and c 2 , which 
can be tuned suitably through Feshbach resonance. For this purpose, we perform a Painleve 
singularity structure analysis to the above fairly generalized system. It is also expected that 
besides BECs the analysis will have wider ramifications in nonlinear optics. 
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This paper is arranged in the following manner. In section II, the three steps involved in 
the Painleve singularity structure analysis, namely the leading order analysis of the Laurent 
expansion in the neighbourhood of a non-characteristic singular manifold, determination of 
the resonances (that is, the powers at which arbitrary functions can occur in the Laurent 
expansion) and analysis of the Laurent expansion for sufficient number of arbitrary functions 
are carried out. It is shown that only for the two specific parametric choices, namely (i) 
c 2 = and (ii) Co = c 2 the system (2) passes the Painleve integrability test. The results are 
analyzed in the final section. 

II. PAINLEVE SINGULARITY STRUCTURE ANALYSIS 

In order to perform the Painleve singularity structure analysis of Eq. (T5]) the dependent 
variables ip±i, an d their complex conjugates are denoted as 



+1 - 



b, ip-i = m, 



n, ip =p, V>o = 9- 



(3) 



Then Eqs. (2) become 



ia t = —a xx + (co + c 2 ) (ab + 2pq)a + (c — c 2 ) mna + 2c 2 np 2 , 

—ib t = -b xx + (co + c 2 ) (ab + 2pq)b + (c - c 2 ) mnb + 2c 2 mq 2 , 

im t = -m xx + (cq + c 2 ) (mn + 2pq)m + (c - c 2 ) abm + 2c 2 bp 2 

—in t = —n xx + (c + c 2 ) (mn + 2pq)n + (c — c 2 ) abn + 2c 2 aq 2 , 

ipt = -Pxx + 2c p 2 q + (c + c 2 ) (ab + mn)p + 2c 2 qam, 



(4a) 
(4b) 
(4c) 
(4d) 
(4e) 
(4f) 



—iq t = -q xx + 2c q 2 p + (c + c 2 ) (ab + mn)q + 2c 2 pbn. 
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The Painleve singularity structure analysis (of an analytic polynomial differential equation) 
is carried out by seeking a generalized Laurent expansion 32 for the dependent variables 

a = 4> a ^ aj(x, t)(fi j , a ^ 0, (5a) 

3=0 

b = ^J2 b ^^^ b ^°> ( 5b ) 

3=0 

m = (jf m,j(x, t)0 J , m 7^ 0, (5c) 
j=o 

n = (f) s ^2n j (x,t)(l> i , n ^0, (5d) 
j=o 

P = <f> € ^2pj(x,tW, Po ^ 0, (5e) 

3=0 

<z = <rJ>0M)^', *>^o, ( 5f ) 

3=0 

in the neighbourhood of the non-characteristic singular manifold (j>(x,t) = 0, with nonvan- 
ishing derivatives <f> x (x,t) ^ and (f> t (x,t) ^ 0. 



A. Leading order analysis 

The leading order behaviour of the solution is analyzed by assuming the forms 

a a (f) a , b ~ & o 0^5 m ~ ^o0 7 ; n ~ n o0 <5 ; P ~ Po4> e ■, 9 ~ 900^ (6) 

for the dependent variables, where a, (3, 7, 5, e and ui are integers to be determined. After 
substituting these forms into Eq. (4) and by balancing the most dominant terms, at the 
leading order one obtains 

a = P = r y = 5 = e = u = — 1, (7) 

with a set of relations 

pi = a m , ql = b n , (8a) 

0x = C ° 2 ° 2 (aAo&o + V m o^o) • (8b) 

Note that there are six functions a , &o> m o, n o, Po an d ?o (besides the arbitrary manifold 
<f>(x,t)) and the above three conditions mean that three of them are arbitrary at this stage 
of the analysis. 
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B. Resonances 



The second step in the singularity structure analysis is to determine the resonances (pow- 
ers) at which arbitrary functions can enter into the Laurent series To obtain the reso- 
nance values, we substitute the following expressions into Eqs. (4) 



a = ao<j) + • • ■ + a 



b = h 



W 1 + ■ ■ ■ + b j( f) J -\ 



m = rriQ(j) 1 + • • ■ + m,-^ 3 1 , n = no(p 1 + ■ • • + rij^ , 



V = Po<P 1 H \-Pj<P 3 \ 



q = qo<P 1 H h <?j<^ 



(9) 



and determine the possible values of j. By collecting the coefficients of J 3 , one can obtain 
a system of six algebraic equations which can be casted as 



D X J 



0. 



(10a) 



where the superscript 'T' denotes the transpose of the matrix and the matrices X and D 
are given by 



X = (aj hj nij n, Pj q s 



(10b) 



D 



( Qi t x al r 2 n a rim a 2r 1 a q 2r 1 p Q a G ^ 

+4c 2 n p 

Tx b 2 Q Q x ri6 n r 2 m b 2r 1 q b 2r 1 p b 

+4c 2 m qo 

r 2 b mo ria m Q 2 nm^ 2rig ^o 2rip ^o 

+4c 2 b p 

rib n r 2 a n nng Q 2 2rxgo^o 2rip ^o 

+4c 2 a g 

nboPa na p ^l^oPo rim Q p <?3 2r 1 a m 
+2c 2 q m +2c 2 a q 
nboqo naoqo rin g rim go 2ri6 ^o <?3 
y +2c 2 p n +2c 2 p b J 



(10c) 
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in which 



n = c + c 2 , r 2 = c - c 2 , 

Qi = -(j - ^ 2 )<^ + 2r i ( fl o&o + Po<?o) + r 2 m n , 

Q2 = -(j - l)(j - 2 )0* + 2r i ( m o«o + Po<?o) + r 2 a b , 

Q3 = -(j - l)(j - 2)(j) 2 x + n (a b + m n ) + 4c p go 

and is a (6 x 1) null matrix. By requiring the determinant of the matrix D to be zero the 
following resonance equation is obtained. 

f(j + l)(j - 3) 3 (j - 4) (4c 2 - 3j(co + c 2 ) + j 2 (c + c 2 )) 2 = 0. (11) 

From Eq. (11) the values of j are obtained as 

j = -1,0,0,0,3,3,3,4,^,^,^,^, (12a) 

where 

1 ( 3(c + c 2 ) + y/9cg + 2c c 2 - 7c 2 2 

■2 \ C0 + C2 " 




(12c) 

All the resonances should be integers for the system (2) to satisfy the Painleve property 
so that movable algebraic branching type critical singular manifolds are avoided. Hence by 
requiring iV\ and iv^ to be integers we find the following two cases: 

Case(i) : c = - ( H 4 — ) c 2 , m = 1, 2, 4, 5, 6, . . . (m ^ 3) 

\ m(m — 3) / 

Case(ii) : c 2 = 0, m = 0, 3. 

Then the integer resonances for both the cases can be written as 

j — — 1, 0, 0, 0, 3, 3, 3, 4, m, m, 3 — m, 3 — m, m = 0,1,2,... (13) 

Note that in the above, j = — 1 corresponds to the arbitrariness of the non-characteristic 
manifold <f>(x, t). 



7 



Case (i): 

In this case for the choices m = 1 and m = 2, we get Co = C2 = c (say), where c is a real 
constant, all the resonances are positive integers (except for j = —1) and are given below. 

j = -1,0,0,0,1,1,2,2,3,3,3,4. (14) 

However for m > 4, the presence of more negative resonances indicates that there may 
exist only particular solutions with lesser number of arbitrary functions in the Laurent ex- 
pansion. For example, the choice m = 4 corresponding to Co = — 2c2, has the resonances 
j = —1,-1,-1, 0, 0, 0, 3, 3, 3, 4, 4, 4. This may be an indication that the Laurent expansion 
(5) does not correspond to a general solution with required number of arbitrary functions, 
but represents only a particular solution. We feel that these choices with m > 4 are the 
candidates for further deeper analysis mathematically on Laurent expansions in the negative 
powers. We do not pursue this problem further here. However, one can perform for example 
a study on the modulation instability of system (2) for these choices of m(> 4) and look for 
solitary wave type solutions which could be of specific physical interest. So, hereafter we 
will consider only the case having resonances (14) with Co = = c. 
Case (ii): 

For the values m = and m = 3, we require C2 = and the system (2) re- 
duces to a set of three coupled nonlinear Schrodinger equations with resonances j = 
— 1,0,0,0,0,0,3,3,3,3,3,4, whose integrability and Painleve analysis have already been 
studied in detail in Ref. 12. So we will not consider this case any further. 

C. Analysis for arbitrary functions 

The next step in the Painleve singularity structure analysis is to show that there exist 
sufficient number of arbitrary functions in the Laurent expansion (5) which can arise at 
the resonance values given by (14) without the introduction of movable critical singular 
manifolds. To prove this, we expand the dependent variables in Eqs. (4) (upto the highest 



8 



resonance value in (14)) as below: 



a = — + a\ + d2<\) + a30 2 + a^cf) 3 , (15a) 


b = ^ + 6i + M> + 6 3 2 + &40 3 , (15b) 

m = — - + mi + m2(p + m 3 (f) 2 + m^ 3 , (15c) 



n = + "i + n 2 + ".30 2 + ^40 3 , (15d) 



p = +Pi +P20 + P30 2 +P40 3 , (15e) 



? = ? + gi + ^0 + g 3 2 + g 4 3 , (I5f) 



where a,j, bj, rrij, rij, pj, qj, j — 0, 1, ... ,4, are functions of (a;, t) to be determined. Then 
by collecting various powers of <ft, we explicitly show that there exist sufficient number of 
arbitrary functions at each index of the resonance values given in (14). As noted above, the 
resonance at j = — 1 corresponds to the arbitrariness of the non-characteristic manifold 0. 



1. Coefficients of (p 3 : 

The set of algebraic equations resulting at this order is 

-2(j) 2 x + 2c(a b + 2p q ) + 2cn pl/a = 0, 

-2<fil + 2c(a 6 + 2p q ) + 2cm ql/b = 0, 

-2(p 2 x + 2c(m n + 2p q ) + 2cb pl/m = 0, 

-2(p 2 x + 2c(m n + 2p q ) + 2ca ql/n = 0, 

-2(fil + 2cpoqo + 2c(a b + m n ) + 2ca q m /p = 0, 

-2(p 2 x + 2cp q + 2c(a b + m n ) + 2cb p n /q = 0. (16a) 

Solving Eqs. (I16al) again results in the already deduced relations (8), 

pi = a m , ql = b n , (16b) 
4>l = c (y a b + V^o^o) • (16c) 

This clearly shows that three out of the six functions (a , bo, m , no, po and go) are 
arbitrary at the triple resonance j = 0, 0, 0. 
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2. Coefficients of 4> 



-2 



At the power 2 , we obtain the following set of algebraic equations expressed in the 
matrix form, 



-i<k) Yi T , 



(17a) 



where 



4c/ 2 

4cp <?o 

4c/! 

4cp go 



4q» <?o \ 
4c/ 2 

4cp <7o 
4c/i 

4cp <7o 



(17b) 



-( 



4c/i 2ca b 2cm n 

2caobo Acli 2cmon 

2ca 6o 4c/ 2 2cm n 
2cao&o 2cm n 4c/ 2 

2c/ 2 2cao&o 2c/ 2 2cm n 2c(/i + Z 2 ) 

y 2ca b 2d\ 2cm n 2cl 2 4cp q 2c(/i + l 2 ) J 

mi ni pi 21 j 
i mo no PO 90 / ' 

( 1 -1 1-1 1 -l) • 

In the above matrix D x we have introduced the quantities l\ and l 2 , which are defined as 

h = a b + p qo and l 2 = m n + p q . (17e) 

In order to make the calculations simpler here and in the subsequent analysis, we use 
the Kruskal ansatz 7 by assuming the singular manifold function <f>(x, t) in the form <f>(x, t) = 
x + p(t), with p an arbitrary analytic function and the aj, bj, rrij, rij, pj, qj are functions of 
t only. One can solve the above six algebraic equations (17) given in the matrix form and 
obtain 



(17c) 
(17d) 



Oi 

a 

b 
mi 
m 
ni 
n 



_Po% 
a b 

a b 
1 



m pi 

2c l 2 po 



ipt 
2c U 



qi 



l 2 \2c po 

L(i£i-i Si 

l 2 \ 2c 1 q 



(18a) 
(18b) 
(18c) 
(18d) 



From equations (18), we observe that the two functions (px and qi) out of the six functions 
ai, bi, mi, ni, pi and q 1 are arbitrary. Naturally, these are associated with the double 
resonance at j = 1,1. 
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3. Coefficients of tp 1 : 
At this order we obtain 

Dx X 2 T = Y 2 T , (19a) 
where the matrix D x being defined in Eq. (117bl) and 



X a = a sa a a fi , (19b) 

\ ao bo mo "0 Po <?o / v ' 

Y 2 = fy« yf yf yf yf yf) ■ (19c) 



Here the elements of Y 2 are given by 

(1) ?:a 0t 2C , 2 2 t \ 

2/2 = \Pi n o + a i b o) , 

ao ao 

4c 

(aiPi9o + aiMo + "iPiPo + ai9iPo + PiSiOo) , (20a) 

a 

(2) z&ot 2c / 2 2 \ 
2/2 = — 7 t (gi^o + Vo) 

»o »o 

4c 

— 7- (&i?iPo + aiMo + "iigigo + OiPigo + Pi?i&o) , (20b) 

(3) _ ^Oi C , 2 + p 2 fe \ 

mo mo 
4c 

(migipo + m x n x mQ + mipig + »iPi£>o + Piqmo) , (20c) 

(4) ^ot 2c / 2 2 \ 
2/2 = [nxuio + (fido) 

n n 

4c 

(ai9i9o + m^rio + nigiPo + PiOin Q + 7iiPig ) , (20d) 

n 

(5) ipot 2c , 2 v 
2/2 = (2pigip + a^mo + m^po + a^po + p-^qo) 

Po Po 

2c 

(aipi&o + OiPia + nipim + aimiq + m^ao + mipira ) , (20e) 

Po 

(6) ^o* 2c / 2 v 
2/2 = (2pigig + nxqimo + m^qo + aiOig + <?iPo ) 

go ?o 

2c 

(bipino + oi^ipo + a^bo + m^no + 6igia + niPi&o) • (20f) 

Proceeding further as in the case of j — 1 and by incorporating the results of j — and 
j = 1 , we express the four functions a 2 , & 2 , wi2 an d ra 2 in terms of the remaining two unknown 
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functions p 2 and q 2 : 



where 



m 2 


n 2 


mo 


no 


n 2 _ 


h 


n 


h 



^2 VPo 9o/ 2c l 2 4 c/i/ 2 

g 2 \ , a 6 (3yf+yf -2^ 1) )+2p go(2yf-yf 



W ' !2cM 2 ' (21d) 

with /i and / 2 being defined in Eq. (17e) and y 2 ^s, j = 1,2, ...,6, are given in Eqs. (20). 
From the above equations (21), we can easily see that two {p 2 and q 2 ) out of the six functions 
a 2 , b 2 , m 2 , n 2 , p 2 and q 2 are arbitrary, as required by the existence of arbitrary functions at 
the double resonance j = 2, 2. 

4- Zeroth order in <f> : 

Collecting now the coefficients at the zeroth order, that is 0°, we obtain 

D 3 X 3 T = Y 3 T , (22a) 
where the matrix D 3 = D x — 21, I is a (6 x 6) identity matrix and 

Xo = ( m M a !M » s ] (22b) 

\ a b m n po qo ) ' 

= (vP yf ] yf ] y? vP yf ] ) ■ ^) 



Y 
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Here 



(1) .fait + a 2 Pt\ J,? . 

y 3 = M (a 2 Pi<?o + aia 2 + a 2 b 1 a + p 2 q\a 

\ a J a 

+n 2 piPo + n x p 2 p Q + a ± q 2 p + axp 2 q + ai& 2 a + aipigi 

2c 

+a 2 giPo + Pig 2 a + Pip 2 n ) {a\bi + n^) , (22d) 

a 

( 2 ) • f b u + hpt\ 4c 

2/3 = ( — ^ ) _ y ( lP2q ° + &2 ^!* + b lP^l + P 2 <?1»0 

+a 2 bib + piq 2 b + aiMo + rn 2 qiq Q + m x q 2 q Q + gig 2 m 

2c 

+&i& 2 a + + b 2 qip ) - — (a^j + miqf) , (22e) 

_( 3 ) _ . ( m lt + m 2 p t \ _ 4c 
m / m 
+m 1 p 1 q 1 + p 2 qim + hp 2 p + b 2 pip + piq 2 m + m x p 2 qo 

2c 

+m 2 p 1 g + mi<7 2 Po + m 2 qip ) + to^tii) , (22f ) 

m 

(4) . ( n u + n 2 p t \ 4c 
y 3 = ~M {n\p 2 qo + n 2 piq + nig 2 p + aig 2 g 

+m 2 nin + n 2 qip + qiq 2 a + qia 2 q + min 2 n + Piq 2 n 



y\ — i ( — J - — (pip 2 b + mim 2 riQ + min 2 m + m 2 nim 



2c 

+P2qin + nin 2 m + 7iiPi<?i) (a ± ql + njmA , (22g) 

n 

(5) .fPit+P2pt\ 4c 2c 

y 3 = M + P2qiPo) {pmqo + m 2 pin Q + minipi 

\ Po J Po Po 
+a 1 b 1 p 1 + a 2 bip + m^po + m 2 nip + ai6 2 p + &i£> 2 ao + ^ip 2 ^o 

+n 2 pim + riip 2 m + 6 2 pia + aip 2 &o + a 2 £>i&o + miq 2 a 

+m 2 gia + aim 2 q + a^m^qx + a 2 mig + a 2 qim + aig 2 m ) , (22h) 

(6) . / git + g2gA 4c \ ^c . 

2/3 = ~M (Pi<?2<7o + p 2 <?i<?o) (<7i<?2Po + m 2 qin + m x n x q\ 

V * / go po 

+ai&igi + 6ig 2 a + miq 2 n + n 2 qim + riiq 2 m + 6 2 gia + min 2 g 
+a 2 fri<?o + ^ 2 nig + aife 2 g + &i^ 2 Po + aig 2 6 + a 2 gi& 

+6 2 nip + n 2 pib + hnipi + riip 2 b + hpxno + 6ip 2 n ) . (22i) 

After a straightforward but lengthy algebra one can solve the above Eqs. (22) and deduce 
the following three expressions: 

'{Acl 2 - 2)g - (4c(/ 2 - h) - 2)g + y« - yf ) 



03 

a (4c/i - 2) 

63 = ((4c/ 2 -2)^-(4c(/ 2 -/ 1 )-2)| + yf -y( 4) ) 
b ~ (4c/! -2) 
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(23a) 
(23b) 



where 



m 3 n 3 



(3) 



) 



+ 



^0(2/3 ~ V3 
2c(h - h) 2 




(23c) 



m n 



2c(h - l 2 ) 



The above expressions (23) indicate the arbitrariness of three functions (n 3 , p 3 and q 3 ) out 
of the six functions a 3 , b 3 , m 3 , n 3 , p 3 and g 3 . Thus system (2) with c = c 2 = c, satisfies the 
requirement of the presence of three arbitrary functions corresponding to the triple resonance 
at j = 3, 3, 3. 

5. Coefficients of (f) 1 : 

In a similar manner, after a lengthy algebra carried out using Maple we have verified that 
the resulting six algebraic equations at the coefficient of 0, which we do not present here for 
want of space, reduce to five equations with six unknown functions. Thus we observe that 
there exists one arbitrary function corresponding to the resonance j = 4, as required. 

Our preceding analysis shows that for the choice Co = c 2 = c, there exist sufficient number 
of arbitrary functions at the resonance values given by Eq. (14). One can proceed further 
to obtain the higher order coefficient functions for all j > 4 in terms of the previous coef- 
ficients without the introduction of any movable critical singular manifold into the Laurent 
expansion for the case c = c 2 = c. So we conclude that the system of three component 
Gross- Pit aevskii (GP) type equation (2) passes the Painleve test only for the two cases (i) 
c 2 = and (ii) cq = c 2 = c and is expected to be integrable. Of course, this fact has 
already been shown for the case c 2 = through the Painleve analysis 12 and the second case 
Co = c 2 = c can be reduced to the (2 x 2) matrix NLS equation 28 which is integrable through 
the 1ST method 33 ' 34 . 

III. CONCLUSION 

In this paper, we have studied the integrability property of the three component Gross- 
Pitaevskii (GP) type equations arising as the evolution equations for spinor condensates by 
applying the Painleve singularity structure analysis. We have identified that only for the 
following two choices of the effective one dimensional coupling constants, (i) c 2 = and (ii) 
Co = c 2 = c, the system (2) passes the Painleve test and possesses Laurent expansion with 
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full complement of arbitrary functions without the introduction of movable critical singular 
manifolds. The integrability of the first choice (c 2 = 0) has been discussed already 12 . For 
the second choice (c = c 2 = c) by applying the 1ST method the multi-bright solitons under 
vanishing 28,29 as well as non-vanishing boundary conditions 31 and multi-dark solitons 30 have 
been obtained by reducing the system (2) to the known 1ST integrable (2 x 2) matrix 
nonlinear Schrodinger equation 33 ' 34 . Our present analysis also shows that the system (2) 
passes the Painleve test for integrability when Co and c 2 are equal and non-zero, in addition 
to the choice c 2 = 0. Apart from finding the choices for which the system (2) can be 
integrable, one can also obtain information regarding the Hirota's bilinearization for such 
cases from the above analysis. This has already been exploited for the case c 2 = to obtain 
multi-soliton solutions 35 ' 36 . Work is now in progress to make a similar analysis for the case 
c = c 2 = c. Also it is of future interest to identify multicomponent integrable systems with 
higher degree of hyperfine spin (F > 1), see for example the coupled evolution equations 
given in 37 , for which also the above type of Painleve singularity structure analysis can be 
carried out. 
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